Holonomic Quantum Computation with Electron Spins in Quantum Dots 
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With the help of the spin-orbit interaction, we propose a scheme to perform holonomic single qubit 
gates on the electron spin confined to a quantum dot. The manipulation is done in the absence 
(or presence) of an applied magnetic field. By adiabatic changing the position of the confinement 
potential, one can rotate the spin state of the electron around the Bloch sphere in semiconductor 
heterostructures. The dynamics of the system is equivalent to employing an effective non-Abelian 
gauge potential whose structure depends on the type of the spin-orbit interaction. As an example, 
we find an analytic expression for the electron spin dynamics when the dot is moved around a 
circular path (with radius R) on the two dimensional electron gas (2DEG), and show that all single 
qubit gates can be realized by tuning the radius and orientation of the circular paths. Moreover, 
using the Heisenberg exchange interaction, we demonstrate how one can generate two-qubit gates 
by bringing two quantum dots near each other, yielding a scalable scheme to perform quantum 
computing on arbitrary N qubits. This proposal shows a way of realizing holonomic quantum 
computers in solid-state systems. 

PACS numbers: 



I. INTRODUCTION 



The emergence of geometrical phases in quantum me- 
chanical systems and their physical and geometrical con- 
sequences were first recognized by Berry and Simon in 
their works on cyclic quantum evolution 0, Soon 
after that, Wilczek and Zee discovered the connection 
between (non-Abelian) gauge fields and the adiabatic 
dynamics of such systems, where the dimension of non- 
Abelian geometric phases is given by the n-fold degener- 
acy of the eigenstates of the Hamiltonian Moreover, 
Aharonov and Anandan Q generalized Berry's idea to 
non-adiabatic evolutions, however, the resulting geomet- 
ric phase is not given anymore by the holonomy in the 
parameter space of the Hamiltonian but in the projective 
Hilbert space. Although, at the time, geometrical phases 
were already known in classical systems [f|, their quan- 
tum mechanical counterparts are physically richer and 
more subtle. 

Based on the above mentioned ideas, a variety of 
schemes for holonomic (HQC) and geometric (GQC) 
quantum computation have been proposed, which re- 
cently attained a considerable attention, and believed to 
be promising candidates to implement quantum comput- 
ers using topological transformations as qubit gates 
In HQC, for example, one can perform quantum com- 
puting by encoding quantum information in the degener- 
ate levels of the Hamiltonian and adiabatically traversing 
closed loops (holonomies) in the parameter space of the 
Hamiltonian. So far, many theoretical (and experimen- 
tal) investigations have been made to implement such 
Hamiltonians in physical systems, for example, confined 
ions in a linear Pauli trap @. The experiment by Jones 
et al. [9] was the first attempt in this direction where 



they realized geometric two-qubit gates between a pair 
of nuclear spins. Note that geometrical phases are gener- 
ally small compared to dynamical phases. Being a small 
effect on top of a large effect, makes it challenging for the 
experimentalists to identify and employ them for quan- 
tum computation. 

Among several proposals for HQC and GQC, solid 
state matrix is usually more desirable due to its potential 
in realizing large scale qubit systems. Specifically, spin 
of an electron in a quantum dot, as a two level system, 
has been shown to be a suitable qubit [Tfj| . meanwhile, 
rapid experimental progress in the field of semiconductor 
spintronics made it possible to access individual electron 
spin in low dimensional systems [ill ]. Manipulating the 
spin of the electrons/holes in semiconductors is, there- 
fore, one of the objectives of spintronics [ll|, tHI- Among 
different tools to achieve this goal is to apply an exter- 
nal magnetic field, in combination with the spin-orbit 
interaction, in a controlled way. Recently, there has been 
a great progress in developing techniques to manipulate 
electrically the electron/hole spins in two dimensional 
electron/hole gases (2DEGs/2DHGs) and quantum dots 
(QDs) |13l - [22l ]. However, most of the previous works on 
confined electrons are based on the assumption that the 
quantum dot itself is almost frozen in real space. More- 
over, the presence of an applied magnetic field is usu- 
ally assumed in order to break the time reversal symme- 
try, which is essential in electron spin resonance (ESR) 
schemes. The question is: what is the dynamics of the 
spin sector if we move the quantum dot in the absence (or 
presence) of the magnetic field? If time reversal symme- 
try is not broken, a convenient way to study the dynamics 
of such a system is to employ non-Abelian gauge fields [3] . 
In spite of the fact that the origin of non-Abelian gauge 
fields in classical/quantum field theories and the current 
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problem is somewhat different, as long as the dynamics 
is concerned, they play the same role. 

Using an effective non-Abelian gauge potential 0, 
[23j . we propose a novel technique to manipulate topo- 
logically the spin of an electron inside a quantum dot 
without using any applied magnetic field. Although the 
rotation of the the electron spin in a moving quantum 
dot has been studied in previous works [12| . |23| , our goal 
here is to implement systematically all necessary single 
qubit gates for quantum information processing. We con- 
sider a setup, where the quantum dot can be moved on 
the substrate at distances comparable to \so, the spin- 
orbit length. In addition, we assume that the electron is 
strongly confined to the quantum dot at a length scale 
Ad <C Aso, where Xd is the dot size . We study the 
case which the confining potential is only displaced par- 
allel to itself by a vector ro (t) without changing its shape 
(see Fig. [2]). Moreover, we assume that the driving elec- 
tric field is classical and only quantize the electron dy- 
namics (for the discussion of the quantum fluctuations of 
the electromagnetic fields which lead to the decoherence 
of the Kramers doublets, see Ref. [23j]). For a moving 
quantum dot around a circular path, we derive an exact 
solution of the Schrodinger equation for the spin sector, 
in the first order in spin-orbit interaction. In addition, 
we show how one can generate different single- and two- 
qubit rotations and perform quantum computing on N 
spins. 



II. THE MODEL AND BASIC RELATIONS 



Quantum Dot 



We consider a lateral quantum dot [24j formed by de- 
pleting the 2DEG via a set of metallic gates that allow 
the quantum dot position ro to be changed at will to dis- 
tances comparable to the spin-orbit length in the 2DEG. 
Such rolling quantum dots can be defined using, for ex- 
ample, a set of gates shown in Fig. [TJ Two layers of 
finger-like gates (separated by an insulator) form a grid, 
which construct the dot confining potential at virtually 
any position under the grid by simultaneously pulsing 
several gates. A relatively different setup is shown in 
Fig. [3j which makes use of a quantum ring and allows 
the quantum dot to be moved along a circular trajectory. 

The electron motion in the plane of the 2DEG, and in 
the presence of a time-dependent dot confining potential, 
is governed by the Hamiltonian 

H(t)=H d (t) + H z + H so , (1) 
where Hd(t) describes the moving dot with one electron, 



H d (t) 



2m. 



+ U(r-r {t)), 



(2) 



with p = —ihd/dr + (e/c)A(r) and r = (x,y) being the 
electron momentum and coordinates, respectively. For 
the vector potential A(r), we choose the cylindric gauge, 
A(r) — B z (—y/2, x/2, 0), where B z is the component of 




FIG. 1: (Color online) A set of metallic gates, deposited on 
top of the heterostructure, control the quantum dot position 
in the (x,y) -plane. The gates that define the current dot po- 
sition are highlighted in a darker color. The superimposed 
gates are separated from each other by an insulating layer. 
The measurement site is used to initialize and read out the 
spin state of the quantum dot with the help of additional con- 
trols (not shown here). By applying time-dependent voltages 
to the gates, one is able to move the quantum dot along a 
desired trajectory, thus forming, for instance, a Wilson loop 
(dashed line). 



the magnetic field normal to the 2DEG plane. The dot 
confinement potential U(r,t) — U (r — ro(t)) changes 
adiabatically with respect to the size-quantization energy 
in the dot and, thus, no transitions between orbital lev- 
els occur. For simplicity, we also assume that the shape 
of the dot confinement does not change in time, while 
the dot is moved along its trajectory ro(t). In Eq. (p}, 
the Zeeman interaction reads Hz — ■ c, with \Ez\ 
being the Zeeman energy and er = (cr x ,a y ,a z ) the Pauli 
matrices. We note that the quantization axis is generally 
not along the magnetic field and one has Ezi = /LtsgyBj, 
where /ib is the Bohr magneton, gy is the g-factor ten- 
sor in the 2DEG, and B is the magnetic field. The last 
term in Eq. ([TJ) , Hso, denotes the spin-orbit interaction 
which has the following general form Hso = h Mp) ' a i 
where h{p) — —h(—p) is an odd-power polynomial in 
p. In GaAs 2DEG with the [001] growth direction, for 
example, the leading order (lowest power in p) spin-orbit 
interaction terms read 

Hso = a(p x a y - p y cr x ) + /3(~p x a x + p y a y ), (3) 

where a and /? are the Rashba and Dresselhaus coupling 
constants, respectively [H, [2(| . 

Considering first a stationary quantum dot potential 
U(r), we set r {t) — > in Eq. (0) and denote the time- 
independent Hamiltonian by H with the following eigen- 
value equation, H\ip ns ) — E ns \ip ns }, where n = 0, 1, 2 . . . 
and s = ±1/2 are the orbital and spin quantum numbers, 
respectively. 

The eigenstates of H are related to the eigenstates of 



Hd + Hz via a unitary transformation (see Appendix) 



(a) 



\4>Ti 



\^n) \Xs 



(4) 



where \ip n ) and \x s ) are obtained by solving Hd\4> n ) — 
E n \tp n ) and H z \Xs) = sE z \Xs), respectively. 

We are mainly interested here in spin dynamics in the 
absence of the applied magnetic field, therefore, we set 
A(r) = 0, and focus on the lowest-in-energy subspace 
n — (for the discussion of the spin dynamics in the 
presence of an applied magnetic field, see Appendix). 
The suitable qubit is then defined as |f) = 1^0,1/2) an d 
\i) = 1^0,-1/2)- In the absence of magnetic fields, the 
quantization axis can, therefore, be chosen arbitrarily. 
However, once it is chosen, all subsequent spin rotations 
are then with respect to this axis. 

Given a linear in momentum spin-orbit interaction in 
Eq. ([3]), and a symmetric confining potential U(r) = 
U(—r), we have explicitly found the generator of the ro- 
tation for the Kramers doublet (see Appendix), as one 
moves the dot along a given path 
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where u are Pauli matrices acting on the Kramers dou- 
blet states and S = icr ■ X^q • r [l6| . According to 
Eq. (JSJ) , the electron state, in the space of a Kramers 
doublet, is rotated during the displacement by an angle 
~ 5ro/Xso- This interpretation of the spin-orbit inter- 
action effect is identical to the standard interpretation 
given to semiclassical electrons. In the latter, the elec- 
tron with the momentum p travels a distance I = pt/m e 
during a time t and changes its spin by an angle propor- 
tional to l/Xso- This coincidence is not accidental, be- 
cause in the semiclassical picture one also assumes that 
the electron moves in a wave packet of an extension that 
is much smaller than Aso- The speed at which the elec- 
tron moves is unimportant at B — 0, because the path 
vo(t) is the only information that determines the spin ro- 
tation. In particular, if the electron travels along some 
path forward and then returns the same way, but not 
necessarily at the same speed, then the initial and final 
spin states coincide. 

We note that, for the linear in momentum Hso, the 
tensor Ag in Eq. ^ is independent of the electron 
orbital state IVVi), which means Eq. (JS|) is valid for all 
symmetric wave packets at the zeroth order of Xd/Xso- 
Therefore, we can consider also a point-like electron, for 
which the orbital wave function reads ip ro (r) — d(r — ro). 
Integrating Eq. ([5]) over an arbitrary path we obtain an 
exact expression for this case, 



so-«-*6(r-r ) X .{<r), (7) 



where the exponent is ordered (to the left) along the path 
of integration. To simplify notations, we use J dr<j to 
denote the contour integral dr' and agree that any 
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FIG. 2: (Color online) (a) Trajectory of the quantum dot cen- 
ter ro(t) and the evolution of the spin state due to displace- 
ment. The spin state changes when going from ro to ro + 5ro, 
due to the infinitesimal transformation (|A.17|I . Since the di- 
rections of Sro can be different in different parts of the curve, 
the infinitesimal transformations do not commute with each 
other and have to be ordered along the path of integration, 
(b) Two oriented circles with the same radii R but different 
centers (with their corresponding holonomic rotations at P) 
show the typical paths for the quantum dot on the 2DEG. (c) 
Moving the dot along the circles will rotate the spin of the 
electron around the 77 axis. For instance, blue (1) and red (2) 
circles in (b) correspond to cf> — and <j> = — respectively. 



exponent of an integral to be ordered along the path of 
integration. The radius vector r (t) gives us the path 
where we choose the beginning of the path at ro = 0, 
and denote the running (present) point of the path by 
ro. Eq. ([7]) determines how the spin of an electron is 
transformed (at B = 0) as the electron is moved along 
an arbitrary path. The difference between the spin and 
the Kramers doublet disappears here, since for point-like 
electrons we can take S — > 0. However, the transfor- 
mation rule in Eq. ([7]) arises from the fact that dS/dr 
remains constant while taking S — > 0. Moreover, it holds 
exactly for the linear in momentum spin-orbit interac- 
tion, because Ago -1 is independent of the orbital state. 
For the p 3 terms, l/Xso is proportional to the electron 
energy and Eq. (J7J can be written only if Hso is first 
linearized around a given energy. 

We developed a code to calculate numerically the spin 
dynamics for an arbitrary path, however as an example, 
we consider a point like quantum dot (strong confining 
potential) moving around a circle with radius R (sec Fig. 
[3]). Using Eq. the dynamical equation for the electron 
spin is then given by 
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(8) 
(9) 



where Xs(<P) — («i>W2) is the spinor, cj> is the angle be- 
tween the starting point vector r (0) and the x axis, 
and Mdr is a Hermitian matrix due to both Dressel- 
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^ Quantum Dot 




FIG. 3: (Color online) Same as in Fig. [T] except that the 
quantum dot is defined using a single layer of metallic gates 
deposited on top of a ring. The ring can be obtained out of a 
heterostructure by means of etching or defined with the help 
of an atomic force microscope using the oxidation technique 
of Refs. EHI. 



haus (A4d) and Rashba (M.r) spin orbit interaction 
(e.g. for the only Rashba interaction we have M. dr = 
A4r = —A4* D ). In addition to trivial solutions for a linear 
path in Eq. ([7]), there are also simple analytical solutions 
for any elliptic or hyperbolic path (x 2 /a 2 ± y 2 /b 2 = 1, 
where a and b are conic section parameters), provided 
a/A + = b/\-. We show here the solution for a <f> = 2ir 
rotation, counterclockwise along the blue (1) circle in 
Fig. [21 for Dresselhaus only where Mdr = M-d an d 
A + = A_ = XR 



Xs(2tt) = exp(--rj • tr)x s (0), 

v = 27r(i-i)(lo,i), 

e A 



(10) 
(11) 



= Vl + 4/A 2 . 



Mathematically speaking, transformation (|TU)) is an el- 
ement of the holonomy group of point P, generated by 
the non-Abelian vector potential 0, Hi| where this non- 
Abelian feature is a direct consequence of the Kramers 
degenerate doublets. Therefore, in the absence of any 
time reversal symmetry breaking interactions (like e.g. 
an applied magnetic field or magnetic impurities), the 
Kramers doublets are robust degenerate states and we 
are able to use the non-Abelian feature of the effective 
gauge potential to manipulate the spin. 



SINGLE-QUBIT ROTATIONS: HADAMARD 
AND PHASE GATES 



III. 



The transformation (1101) is nothing but a rotation 
along the vector rj at point P which is in the x-z plane 
and makes an angle 9 = arctan(j) with the z axis, see 
Fig. [2J Therefore, for large values of A (small circles), 



9 « and one would able to realize arbitrary rotations 
around z axis (phase gate) with reasonable precision. 
Moreover, moving counterclockwise along the red (2) cir- 
cle in Fig. [2] leads to the same result as in the blue (2) 
circle, but now the rotation takes place in the y-z plane. 
Therefore, depending on the orientation of the circles and 
their corresponding radii, we can, in principal, achieve all 
kinds of rotations around the Bloch sphere. 

To be more specific, we show how to generate the 
Hadamard gate by rotations around two non-orthogonal 
axes (in our case z and r] directions) (30j and, for conve- 
nience, we only consider the Dresselhaus term (a = 0). 
As shown in Eq. (|10[) . circles with different radii and/or 
orientations lead to different rotations. In particular, if 
we go counterclockwise along a full circle from point P 
which makes an angle <j> with x axis (see Fig. [2]), the 
electron spin will transform as follows 



— cos ire 



■ sm7re, 



Ae 
U22, 



c 

sin 7re, 



U12 = 



where Uij is the unitary transformation which acts on 
the initial spin state. Geometrically, U is the matrix 
corresponding to a rotation around rj axis which lies in 
the same plane as z axis and P (see Fig. [21c). 

Hadamard gate can be achieved (up to a global phase) 
by a clockwise \ rotation around y axis followed by a 
counterclockwise 7r rotation around z axis, 



H 



1 

V2 



1 1 

1 -1 



iU z (n) U y (- 



(12) 



The rotation around z with an arbitrary angle has al- 
ready been discussed above, therefore, we show here how 
one can implement U y (Q). The main problem is that, 
according to Eqs. [1011 1] . the magnitude and the direc- 
tion of the vector rj are not independent variables (both 
are functions of the variable A). Moreover, the vector 
77 is not in general orthogonal to z axis. To obtain an 
arbitrary rotation around, say, y axis, we need to per- 
form 3 rotations (two around rj and one around z axis). 
Specifically, we want to know for which values of A, ±-| 
rotations around y axis can be achieved. One can show 
that [H 



^(7) = U V (9)UMU V (8), 



(13) 



cos rj 



cos 2 9 sin 2 \ ± cos \ J 1 - cot 2 8 sin 2 3r 



cos 2 9 sin 2 ^ — 1 



tan (o) = — — 
l sin 



sin r\ cos f 



cos rj cos z 



V = \m, 



where for our purpose, U y (^), we need to evaluate 
Eq. (fT3")l at 7 = 7r/2. Obviously, there are an infinite 
number of solutions corresponding to different values of 
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FIG. 4: (Color online) Schematic of the two-qubit rotation 
setup. The confined electrons are brought together at the 
intersection of the circular paths (with their potential pro- 
file shown in the inset at this point), where they interact 
via Heisenberg exchange interaction. The (holonomic) single- 
qubit operations are done at position P (P') while the (ex- 
change) two-qubit gates are performed at position G (G'). 



A. Therefore, the tt/2 rotation around y axis, and conse- 
quently the Hadamard gate, is achievable in our scheme. 
Together with the phase gate (arbitrary rotation around 
z axis), all single qubit operations can be realized. 

We note that the quantum fluctuations of the driving 
field will lead to the decoherence of the Kramers doublets 
even in the absence of the applied magnetic field, how- 
ever, this rate saturates at the zero field strength [2J|. For 
typical lateral GaAs quantum dots with the size A^ ~ 50 
nm which corresponds to the orbital quantization ujq ~ 1 
meV, the estimated decoherence rate T ~ fis^ 1 [23[. Ob- 
viously, to implement efficient single qubit gates by using 
a ring of tunnel-coupled quantum dots, the pulsing and 
the total travel time of the quantum dot should be then 
much smaller than the decoherence time On the 

other hand, the adiabaticity criterion puts a limit on the 
velocity v of the moving quantum dot in order to keep 
the electron in its ground state doublet. Therefore, one 
needs to satisfy TXso <C |u| <C Xd^o at an y moment in 
time for a quantum dot which is displaced on the scale 
of Aso • Recent experiments on tunnel-coupled quantum 
dots show the ability to transfer the electron wave func- 
tion over few hundreds of nm (from one dot to the other) 
in 1 ns |2i| . We observe that, for instance, for the spin- 
orbit length Aso ~ 3 /im, this pulsing time is perfectly 
within the range of the above mentioned condition and 
our scheme is, therefore, experimentally feasible. 



from the position P (P') to G (£?'), as shown in Fig. [4] 
The top gates are introduced to control over the wave 
functions of the confined electrons at the touching points 
of the circles. By lowering the potential barrier (V g ) be- 
tween two quantum quantum dots, the residing electron 
spins can couple to each other (due to the overlap of their 
orbital wave functions) via the Heisenberg exchange in- 
teraction [1(1 Hlj]. It nas already been shown that by 
electrically engineering the gate potentials, one can gen- 
erate the SWAP gate, and eventually the CNOT gate, 
between two spins [1(|. However, this additional step, 
i.e. moving two quantum dots towards each other to 
perform two-qubit gates, leads to a spin rotation of each 
electron, see Eqs. (|8l9p , and therefore this partial spin 
dynamics should also be taken into account. For practical 
purposes, we assume that the radius of the (holonomic) 
single-qubit gates (R) is smaller than the radius of the 
two-qubit circles (R ! ), to avoid spatial overlap of different 
local single-qubit operations, see Figs. [Ubgl The read- 
out part of the scheme is accomplished by applying an 
external magnetic field, using different techniques such 
as spin-to-charge conversion [24j |. 

As an example, we observe that our scheme is able 
to produce the so called cluster states on N qubits [32| . 
The peculiar properties of cluster states made them a 
suitable candidate for realizing quantum computers in 
quantum optical and solid state systems [32-35]. One of 
the main advantages of one-way quantum computation 
is that these set of entangled states (cluster states) are 
produced once, and then quantum computation is done 
by local (adaptive) measurements of the qubits. There- 
fore, there is no need to perform two-qubit gates during 
the information processing. In Fig. 01 we showed a one 
dimensional linear chain of qubits, however, this scheme 
can be easily generalized to higher dimensional lattices, 
and in principal be used to generate cluster states and 
perform holonomic one-way quantum computing in solid 
state environments. 
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IV. TWO-QUBIT GATES, READ OUT, AND 
QUANTUM COMPUTATION 

For quantum computing, in addition to single qubit 
rotations, two qubit operations are needed as well (e.g. 
CNOT gate). Here, we propose a scalable setup to per- 
form quantum computation on arbitrary N qubits (see. 
Fig. 3]) . To this end, we move two quantum dots around 
a circular path and bring them close to each other, e.g. 



Appendix: The Non-Abelian Spin Rotation 
Generator 

The transformation matrix S = — S' in Eq. ^ can be 
evaluated by perturbation theory in Hso or by diagonal- 
izing the Hamiltonian H for a specific potential U(r). At 
the leading order in H$o, S satisfies the operator equa- 
tion 

[H d + H z ,S]=H so , (A.l) 
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whereas the energy levels of H coincide, at this order, 
with the energy levels of H d + H z . A detailed analysis 
of the transformation in Eq. (|4| is given in Refs. (l6Ll3(ij. 
Next we consider the Schrodingcr equation, 



ih 



d^(r,t) 
dt 



H(t)*(r,t), 



(A.2) 



in the presence of a time-dependent displacement-vector 
ro(t). At each moment in time, the Hamiltonian H(t) 
has an instantaneous basis of states, which we denote 
by l^nsro); where the index ro indicates that the dot 



is centered at r . Obviously, the states |$ r 



can be 



obtained from Eq. Q by means of a displacement by 
the vector — ro. In the presence of a magnetic field, the 
instantaneous eigenstates read 

$ nsro (r) = e^ hc ^^T_ ro ^ ns (r) , (A.3) 

where / (r, ro) is a gauge function satisfying the equation 

df (r,r ) 



dr 



= A{r - r ) - A(r) 



(A.4) 



In Eq. (|A.3|) . T a — exp(a ■ d/dr) denotes the translation 
operator by a vector a and for the cylindric gauge A(r) — 
B z (—y/2, x/2, 0), we can choose / (r, r ) = r ■ A(r). 

The solutions of the Schrodingcr equation (|A.2|) can 
be looked for in terms of the instantaneous basis in Eq. 



tf(r,i) 
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(r,t), (A.5) 



where the coefficients a ns (t) satisfy the normalization 
condition ^2 ns \a ns (t)\ 2 = 1 and we have used the no- 
tation $ nsro ( t ) (r,t) = exp (-i Enst /ti) $ nsro ( t) (r) for 
the Schrodinger picture of Eq. (|A.3j) . By substituting 
Eq. (TOt into Eq. (|A~2l) we obtain a set of equations for 
a ns {t) 



da 



ns 1> 

dT h 



!^o(i)-(^|p(i)|Vv s <)e 1, 



(A, 



where Vo(t) = droit) /dt is the velocity of the slipping 
dot and linen's' = (E ns -E n ' s ') /h. The quantity p(t) 
depends on t only through T*o(i) and is defined as follows 

(A.7) 



P = -inl- e - Tro d M^T_ 

or c 



dr 



For our choice of gauge, i.e. cylindrical gauge, we obtain 
p = —ihd/dr — (e/c)A(r + ro). Note that the choice of 
/(r,ro) is not unique; In general, one can also include 
terms of the form (?o( r o) + s 53( r o) IVVis) (ipns\ and, at B = 
0, additionally terms of the form gi{r ) \ipns)(tpn,-s\ + 
isg2{ro)\^ns)(i'n-s\, where g 3 (r ) are arbitrary real 
functions of ro- 

Next we consider a specific situation for which we can 
further simplify Eq. (|A.6[) . 



We can further define a resting qubit at a position ro 
using the transformation in Eq. (|A.3|) . Let the quantum 
dot be driven along a trajectory r (t) between two points 
r a = and rs during a time interval T, such that 



ro(0) = ta, 
MO) = 0, 



r (T) = r B , 
v o (T) = 0. 



(A.8) 
(A.9) 



The probability for the qubit to leak out of its subspace 
by the end of the pulse is given by 



flcak = \a ns (T)\' 



(A.10) 



s = ±l/2 



The coefficients a ns (T) can, therefore, be found by solv- 
ing Eq. (|A.6|) with the initial condition ^ s |a Os (0)| 2 = 1. 
At the leading order in the driving, we have 



a ns (T) ~ - J dt 



Vo(t) ■ {lpna\p(t) W)e 



, (A.ll) 



where |u) = \ipoa) denotes the qubit state at t = 0. Since 
the matrix elements (ipns \ P{t) |f ) do not depend on time 
for n 7^ 0, the coefficients a ns (T) in Eq. (jA.111) are, thus, 
proportional to the Fourier transform of the quantum 
dot velocity Vq (t) evaluated at the orbital transition fre- 
quency uinsOa — (jJo- It is, therefore, sufficient to devise 
pulses of vo(t) with the spectral weight below the orbital 
frequency ujq in order to avoid leakage from the qubit 
subspace. 

It is convenient to have an adiabaticity criterion based 
on the differential properties of v (t). We note that 
Eq. (|A.6|) can be rewritten in terms of the new unknowns 
a ns = a ns exp (—itE ns /h) as follows 



da 7i 



dt 



,(t)d r< 



(A.12) 



rt r , 



>(t) = E ns 5 ns ,n's> ~ V (t) ■ (lp ns \p{t) \lp n ' s ') . (A.13) 



One can identify Eq. dATTB]) with H(t) = H(t) - ihd/dt, 
expressed in the time-depended basis (|A.3[) . For the case 
of B = 0, Eq. (IA.13|) has been previously obtained in 
Ref. [H|]. The virtue of Eq. (|AT3l is that U nsn > s '{t) 
depend on time only through a perturbation oc Vo(t), 
which vanishes at t — 0,T. Applying the adiabaticity 
criterion to the Hamiltonian in Eq. (|A.13|) for the orbital 
transitions out of the qubit subspace, we obtain that the 
condition 



dv 



p\a) 



(A.14) 



must be satisfied at any moment in time in order for the 
pulse to be adiabatic. 

If r (t) changes adiabatically with respect also to the 
Zeeman energy Ez, then |a s (f)| is independent of time, 
i.e. the qubit follows adiabatically the change of its basis 
states. In the opposite case, when B = 0, the states \i[> na ) 
in Eq. @ are degenerate with respect to the spin index 
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s to all orders of Hso, due to the Kramers theorem. In 
this case, the change of the instantaneous basis can be 
interpreted as a unitary operation on the qubit. In order 
to tell what is the qubit instantaneous basis at B = 0, 
one has, in principle, to consider a finite B and follow the 
energy levels of the quantum dot in the limit of B — > 0. 
Here, it is important to note that the spin-orbit inter- 
action gives rise to an anisotropic Zeeman interaction at 
the second order of Hso [3- As a result, the spin quan- 
tization axis and the magnetic field are not necessarily 
aligned with each other. To avoid the need of state find- 
ing, we denote (ns\e~ 
a n * = 0(H 2 so ). 



Eqs. (|A.18I) and (|A.19[) account also for the Aharonov- 
Bohm phase due to the orbital magnetic field. It is, 
therefore, convenient to subdivide Q into Abelian and 
non-Abelian parts, 



V ss' y.ss' 



P 



(A.20) 
(A.21) 



For a point-like quantum dot, we obtain the Abelian gen- 
^ erators Q a ss , — (—ie/hc)S ss >A(r ), recovering, thus, the 

by a ns and remark that usua l expression for the Aharonov-Bohm phase (e* VAB ) 



Returning now to Eq. (|A.5|) . we consider an infinitesi- 
mal displacement of the quantum dot in the (x,?/)-plane 
by 5ro and derive the corresponding generators of the 
qubit transformation under translations. We encode the 
qubit into the instantaneous states of the n-th orbital 
level of the Hamiltonian ([1]) . 

Let ro(t) — To be the position of the quantum dot 
center at time t and r (t + St) = r + Sr be the new 
position at time t + St. The infinitesimal transformation 
that takes the state 'S(t) to a new state ^>(t + St) is given 
in Eq. (IA.2I) . Starting from a basis state § n s(t) at time 
t, we obtain the following state at time t + St, 



V ns (t + St) = $ ns {t) - -H(t)$ ns (t)5t. 



(A.15) 



The overlap of this state with the basis state <& ns / (t + St), 
generates an infinitesimal transformation of the Kramers 
doublet, where from Eq. (|A.3|) . we find the basis state at 
time t + St, 



$ ns >(t + 5t) = ^ ns '(t) 



dr 
dt 



p + H(t) 



$ns>(t)6t. 

(A.16) 

Thus, the desired infinitesimal transformation reads 

($ ns >{t+St)\V ns (t+5t)) = 5 s , s +^5ro-{$ ns >(t)\p\$ns(t)). 

(A.17) 

For a qubit that is encoded into the instantaneous states 
of the n-th orbital level of the Hamiltonian (p}, the in- 
finitesimal transformation (|A.17j) can be rewritten as 



\s(t)) ->exp(ff.*r )Kt)). 



(A.18) 



Here, \s(t)) denotes the qubit state at time t and the 2x2 
matrices 



^WwIpIVw) 



(A.I9) 



are the corresponding generators of the transformations 
that take place on the qubit under parallel translations of 
the quantum dot on the substrate. In deriving Eq. (|A.19|) 
we made use of our choice of gauge, see the text below 
Eq. CO}- 

It is important to note that, along with spin-orbit 
interaction-induced SU(2) transformations on the qubit, 



</?AB 



■ dr Q 



(A.22) 



in going around a closed path C. Note that it is always 
possible to sum up independently the phase due to the 
Abelian generators, because [Q a ,Q na ] = 0. In what fol- 
lows, we focus on the non-Abelian generators Q na since 
they give rise to useful unitary operations on the qubit. 

To calculate the matrix elements (& ns r(t)\p\Q ns (t)) , 
we make use of Eq. (|A.3|) and the following property 
e -'e/ /ficpgie/ /he _ p an( j obtain that 

($„ s '(*)|p|$« s (f)) = e^ B -'-^>'(V^T ro pT_J^„ s ), 

(A.23) 

where \ipns) are the states in Eq. ((4]) and E ns are the ener- 
gies corresponding to these states. Obviously, if the Zee- 
man energy is large, the exponential factor in Eq. (| A.23|) 
oscillates rapidly as a function of time and the transfor- 
mation in Eq. (IA.17|) averages out to unity. For the latter 
to take place, it is sufficient that 



dr (ip ns \T ra pT_ ro \i/; n _ 6 
dt 



E,, 



En 



< 1. 



(A.24) 



Estimating further \{ip ns \T ro pT_ ro \ip nj _ s )\ ~ h/X S o and 
E ns — E n ,- S 1=3 Ez, we obtain that the spin rotator is 
inefficient at small speeds of the dot, hfo <C EzXso- 

In the absence of magnetic fields, the transformation 
in Eq. (|A.17|) acquires the form 



d 

A = 1 + Sr ■ (Vw k-Kw- 
or 



(A.25) 



Note that Eq. (|A.25|) can as well be derived from 
the infinitesimal version of the identity \Q ns (t)) = 
T S r B T-Sr \<S> ns (t)) = T Sro \<$> ns {t + 5t)). Thus, the spin 
rotation takes place (at least at B — 0), because the con- 
finement defines local Kramers states, which differ from 
each other along the dot trajectory. An illustration of the 
dot trajectory is given in Fig. [2] The radius- vector ro(t) 
describes a curve as a function of time and, as the dot is 
moved along that curve, the local Kramers state changes. 
The infinitesimal transformation in Eq. (|A.25[) [or more 
generally in Eq. (|A.17[) ] has to be ordered along the path 
of integration when integrated over Sr . This ordering 
occurs because the spin matrices in Eq. (|A.25[) do not al- 
ways commute with each other at different points of the 
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path due to generally different directions of 5ro at these 
points. 

Further, it is convenient to refer to the Kramers dou- 
blets \$> ns (t)) as to spin states that are locally defined at 
each point of the dot trajectory. Mathematically, we per- 
form a mapping given by the following canonical trans- 
formation, 

\<f> ns ,r }=e ie f/ hc T_ ro e- s \i )n )\ Xs ), (A.26) 

which is obtained by substituting Eq. ((4]) into Eq. (|A.3|) 
and omitting the free evolution factor e~W i > Enst from 
^na{i'i't)- For a given quantum number n, this transfor- 
mation, obviously, maps the Kramers doublet at position 
r onto a spin-1/2 space: |x s ), (s — ±1/2). Equation 
(|A.25[) can then be rewritten in an operator form 

A = l + ( 5r .(^|e s |-e- s |V n ), (A.27) 



where S contains Pauli matrices, which should now be 
regarded as effective operators in the Hilbert space of a 
local Kramers doublet (n, ro). 

For a small quantum dot the transformation matrix S 
is small, because <C \so [36]. In this case, one can 
expand the transformation to the first order, e ±s ~ 1±S. 
Then, Eq. (|A.27[) acquires the form (B = 0) 



A = l-<Sr -(^„| — \ip n ), (A.28) 



where we used (tp n \d / dr\tp n ) = 0. 



M. Berry, Proc. Roy. Aca. A 392 45 (1984). 

B. Simon, Phys. Rev. Lett. 51, 2167 (1983). [19 
F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984). [20 
Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 [21 
(1987). 

A. Shapere and F. Wilczek (eds.), Geometric Phases in [22 
Physics (World Scientific, Singapore, 2002). 
P. Zanardi and M. Rasetti, Phys. Lett. A 264, 94 (1999). [23 
E. Sjoqvist, APS Physics 1, 35 (2008). 
L.-M. Duan, J.I. Cirac, and P. Zoller, Science 292, 1695 [24 
(2001). [25 
J. A. Jones, V. Vedral, A. Ekert, and G. Castagnoli, Na- [26 
ture 403, 869 (2000). [27 
D. Loss and D.P. DiVincenzo, Phys. Rev. A 57, 120 [28 
(1998). [29 

D. D. Awschalom, D. Loss, and N. Samarth (eds.), 
Semiconductor Spintronics and Quantum Computation [30 
(Springer, Berlin, 2002). 

W.A. Coish, V.N. Golovach, J.C. Egues, and D. Loss, [31 
Physica Status Solidi (b) 243, 3658 (2006). 

E. I. Rashba and ALL. Efros, Phys. Rev. Lett. 91, 126405 [32 
(2003). 

E.I Rashba and ALL. Efros, Appl. Phys. Lett. 83, 5295 [33 
(2003). 

M. Duckheim and D. Loss, Nature Physics 2, 195 (2006). [34 
V.N. Golovach, M. Borhani, and D. Loss, Phys. Rev. B [35 
74, 165319 (2006). [36 

C. Flindt, A.S. Sorensen, and K. Flensberg, Phys. Rev. 
Lett. 97, 240501 (2006). 

J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 



Rev. Lett. 78, 1335 (1997). 

G. Salis, et al. Nature 414, 619 (2001). 

Y. Kato, et al. Science 299, 1201 (2003). 

Y. Kato, R.C Myers, A.C. Gossard, and D.D. 

Awschalom, Nature 427, 50 (2004). 

D. V. Bulaev and D. Loss, Phys. Rev. Lett. 98, 097202 
(2007). 

P. San- Jose, B. Scharfenberger, G. Schon, A. Shnirman, 
and G. Zarand, Phys. Rev. B 77, 045305 (2008). 
R. Hanson, et al. Rev. Mod. Phys. 79, 1217 (2007). 

E. I. Rashba, Sov. Phys. Solid State 2, 1109 (1960); 
G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

R. Held, et al. Appl. Phys. Lett. 73, 262 (1998). 

A. Fuhrer, et al. Nature 413, 822 (2001). 

M. Nakahara, Geometry, Topology and Physics (IOP 

Publishing Ltd., 1990). 

R. Hanson and G. Burkard, Phys. Rev. Lett. 98, 050502 
(2007). 

G. Burkard, D. Loss, and D.P. DiVincenzo, Phys. Rev. 
B 59, 2070 (1999). 

H. -J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001). 

R. Raussendorf and H.-J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

M. Borhani and D. Loss, Phys. Rev. A 71, 034308 (2005). 
H.-J Briegel, et al. Nature Physics 5, 19 (2009). 
M. Borhani, V.N. Golovach, and D. Loss, Phys. Rev. B 
73, 155311 (2006). 



